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Fluctuations in nuclear collisions can be measured as a function of momentum-space binning 
scale over a scale interval bounded by detector two-track resolution and acceptance. Fluctuation 
scale dependence is related to two-particle correlations by a Fredholm integral equation. That 
equation can be inverted by standard numerical methods to yield an autocorrelation distribution 
on difference variables as a projection of the full two-particle distribution which retains most of the 
correlation information in a more compact form. Autocorrelation distributions are typically more 
easily interpreted in terms of physical mechanisms than fluctuation measurements. 
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I. INTRODUCTION 

The dynamics of relativistic nuclear collisions are stud- 
ied via the momentum structure of the hadronic fi- 
nal state pj. The correlation structure of pair density 
p(p\,y-i) on two-particle momentum space provides qual- 
itatively new information about collision dynamics com- 
pared to single-particle density p(p). Since the full two- 
particle momentum space cannot be visualized, it is nec- 
essary to project the pair distribution onto subspaces, 
retaining as much of the correlation structure as possible 
in the projections. For relativistic nuclear collisions the 
autocorrelation distribution Q represents a nearly loss- 
less projection Q of the two-particle momentum distri- 
bution onto visualizable subspaces (i.e., projections onto 
two dimensions). 

Measurements of the fluctuations of event-wise multi- 
plicity n and mean transverse momentum (p t ) have been 
proposed to study the QGP phase boundary as probed by 
heavy ion collisions 0, 0, IE 13 • Multiplicity and mean- 
Pt fluctuations may also result from hard parton scatter- 
ing and the hadronization process @, . Direct physical 
interpretation of fluctuation measurements is typically 
difficult. However, the bin-size or scale dependence of 
fluctuations is determined by the correlation structure 
of the underlying distribution 0, ^| . We have deter- 
mined that fluctuation scale dependence can be inverted 
by a numerical procedure to obtain autocorrelation pro- 
jections which can be interpreted directly. 

In this paper we derive an exact relation between 
fluctuation scale dependence and autocorrelation distri- 
butions in the form of an integral equation. We de- 
scribe the procedure required to invert the integral equa- 
tion 01 an d thereby obtain autocorrelations from fluc- 
tuation scale dependence. We apply that procedure to 
Monte Carlo simulations and examine the relation of the 
revealed correlations to known Monte Carlo physics and 
the precision of the autocorrelations obtained by inver- 
sion compared to direct pair analysis. 



II. PREVIOUS WORK 

A. Scaled factorial moments 

Intermittency analysis introduced the concept that 
fluctuation dependence on bin size or scale should be 
examined over a range of scales to study particle pro- 
duction mechanisms in nuclear collisions |l3j| . The prin- 
cipal motivation was the conjecture that some particle 
production mechanisms may produce fluctuations which 
saturate at a particular scale, whereas other mechanisms 
(e.g., random cascades) may exhibit power-law behavior 
limited only by finite experimental resolution. The pro- 
posal confronted the problem of detecting fluctuation ex- 
cess (beyond that due to finite particle number) by intro- 
ducing normalized factorial moments which take on unit 
value for Poisson fluctuations. Power-law (fractal) behav- 
ior is expected theoretically from both QCD branching 
processes and some models of the QCD phase transition 
(e.g., critical fluctuations), typically motivating log-log 
plots of normalized factorial moments vs binning scale. A 
major experimental program to detect and study power- 
law behavior based on factorial moment analysis has met 
with considerable success 111. 



B. Correlation integrals 

The correlation integral [l4|, introduced as a supple- 
ment to factorial moment analysis, establishes the rela- 
tion between moments (ordinary moments, factorial mo- 
ments, cumulants) of a multiparticle density distribution 
and integrals of that distribution over certain momen- 
tum subspaces. The integrand of a correlation integral 
can be expressed as a projection of the full multiparticle 
space onto a subspace of momentum difference variables. 
If that projection is obtained by averaging it is an au- 
tocorrelation. A connection is thus established between 
running integrals of autocorrelations on difference vari- 
ables and scale dependence of distribution moments. 
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C. Heavy ions and the integral equation 

Normalized factorial moments plotted on a log format 
vs logarithmic scale variables have been broadly applied 
to elementary collisions. In a heavy ion context the anal- 
ysis problem changes significantly. The accessible scale 
interval is typically small - 1-2 decades, but the statisti- 
cal power of the data is typically very large and the scale 
dependence may be highly structured over small scale in- 
tervals (as opposed to a fixed power-law trend over a large 
interval). The ensemble-averaged distribution providing 
the reference may be complex, and strongly dependent on 
collision parameters such as energy and centrality. Thus, 
the advantage of normalized factorial moments in terms 
of the Poisson reference does not offset difficulties with 
low- multiplicity data (empty-bin effect). 

Those features lead in heavy ion analysis to three mod- 
ifications of the factorial- moment analysis: 1) ordinary 
moments replace factorial moments, 2) the differential 
fluctuation measure is a linear combination of moments 
(e.g., variance or higher cumulant differences) with the 
cumulant reference provided by local bin averages (as for 
a 'vertical' factorial moment analysis) and 3) the scale de- 
pendence is plotted on a linear format - e.g., linear vari- 
ance difference vs linear bin size. That linear framework 
combined with the correlation integral concept leads to 
the definition of an invertible integral equation connect- 
ing fluctuation scale dependence to autocorrelation dis- 
tributions, as described below. 



III. STATEMENT OF THE PROBLEM 

A. Motivation 

Motivations for this work are three- fold: 1) Fluctua- 
tion measurements have been promoted to study branch- 
ing processes (intermittency, QCD bremsstrahlung) and 
critical charge, flavor and isospin fluctuations associ- 
ated with the QCD phase boundary. Interpretation of 
fluctuation measurements is indirect and model depen- 
dent, whereas inversion to correlations usually leads to 
straightforward interpretations; 2) Direct calculation of 
two-particle correlations is computer intensive - 0(n 2 ), 
whereas calculation of fluctuations is fast - 0(n). If cor- 
relations are derived from fluctuations considerable time 
savings are possible, most important for heavy ion data 
samples of millions of events with thousands of particles; 
3) Two-particle correlations computed directly as pair 
ratios or differences contain an offset ambiguity related 
to relative normalization of object and reference distri- 
butions. The ambiguity is eliminated for correlations ob- 
tained by fluctuation inversion. 



B. Basic concepts 

The single-particle number density p(p) <-> d 3 n/dp 3 
with integral (n) is estimated by accumulated particle 
counts on a binned momentum space (histogram). His- 
tograms for individual events sample the (unobserved) 
event-wise parent density distributions determined by 
collision dynamics. If the parent is the same for all 
events the event-wise estimators deviate from the en- 
semble mean only according to Poisson statistics. The 
ensemble-averaged histogram then contains all informa- 
tion about the parent distribution and the dynamics. 

If the parent distribution for individual events deviates 
significantly from the mean, the variation produces non- 
Poisson fluctuations of histogram bin contents. The scale 
(bin size) for each momentum component determines the 
amount of information in the fluctuations. In this work 
fluctuations are measured by variances computed for bins 
over a range of scales. The scale dependence of variances 
reveals the spatial structure of event-wise changes of the 
parent distribution (as outlined in [13j for factorial mo- 
ments), and therefore the collision dynamics. 

The two-particle pair-number density p(p±,p2) «-> 
d 6 n/dpidp2 with integral (n(n — 1)) is similarly estimated 
by pair counts in a histogram on the two-particle mo- 
mentum space. Just as a single-particle histogram has 
a mean-value component and possibly non-Poisson bin 
fluctuations about the mean, the two-particle distribu- 
tion has similar components. However, the mean values 
of the two-particle histogram are related to the fluctu- 
ations of the single-particle distribution, suggesting an 
algebraic relationship between two-particle correlations 
and single-particle fluctuations. 

C. Correlations and autocorrelations 

Two-particle correlations are measured by the differ- 
ence between object and reference histograms. The ob- 
ject distribution is the ensemble mean of event-wise two- 
particle histograms (sibling pairs). The reference distri- 
bution is either the cartesian product of single-particle 
distributions or a histogram of particle pairs formed from 
different events (mixed pairs). In either case the refer- 
ence represents particles produced as independent sam- 
ples from the single-particle parent. The object con- 
tains additional information on correlated particle pro- 
duction, for instance emission from several relatively- 
moving sources. 

Two-particle difference histograms on momentum 
components pseudorapidity n and azimuth (f> are observed 
to be approximately invariant along the sum axis (e.g., 
Vi + V2) for relativistic nuclear collisions. Such invari- 
ance is consistent with two possibilities: 1) the correla- 
tion mechanism is truly invariant on one or both vari- 
ables (e.g., hydrodynamic flow) or 2) structure is ran- 
domly distributed in the single-particle momentum space 
from event to event, but produces a structure invariant 
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on the sum variable in the ensemble-averaged histogram 
(e.g., jet correlations). In either case, the structure on 
the difference variables (e.g., rj\ — 772) then represents all 
correlation information as a projection by averaging: an 
autocorrelation. 

An autocorrelation compares a distribution f(x) to it- 
self in the form of an average of product f (x\) f (x 2 ) along 
the sum variable x\ + X2 for each value of the differ- 
ence variable X\ — x 2 . That projection by averaging onto 
the difference variable is a lossless projection strategy if 
correlations are invariant on the sum variable (station- 
ary) . Autocorrelations are well-established in time-series 
analysis (where the difference variable is the lag r) [lflj . 
They also appear in various aspects of nuclear collision 
analysis. Direct determination of joint autocorrelations 
on momentum subspace (77, (f>) has been extensively im- 
plemented in multiparticle correlation analysis of fixed- 
target p-p experiments [Tril ]. Autocorrelations also pro- 
vide unbiased access to jet correlations, as we demon- 
strate in the section on Monte Carlo simulations. The 
examples above involve direct determination of autocor- 
relations by pair counting in histograms. Alternatively, 
autocorrelations can be obtained by inverting the bin-size 
or scale dependence of fluctuations. 



D. A simple example 

An event-wise parent distribution on x is defined as 
a gaussian superposed on a flat background. In each 
'event' the gaussian appears at a random position. The 
parent for each event is sampled randomly. The ensemble 
mean sample distribution on x (single-point distribution) 
is uniform, retaining no information about the gaussian 
shape. The sample-pair distribution for each event is 
histogrammed on (a; 1,0:2) (two-point distribution), and 
the object distribution is the ensemble-average pair dis- 
tribution, which includes a uniform component from the 
flat background, a second uniform component from gaus- 
sian pairs distributed randomly on (xi,x 2 ), and a ridge 
along the sum diagonal representing sample pairs from 
the same gaussian (gaussian self pairs - sample self pairs 
are not included in the histogram). The reference dis- 
tribution formed with sample pairs taken from different 
events is uniform (cartesian product of single-point dis- 
tributions). The object-reference difference reveals the 
ridge along the sum diagonal, which represents the gaus- 
sian component of the parent. The projection of that 
difference histogram onto difference variable x\ — x 2 by 
averaging along the sum variable is an autocorrelation. 
The principal feature of the autocorrelation is a gaus- 
sian, with amplitude and width determined by the mean 
of the parent gaussians. 

From the same event ensemble the sample variances 
from a fluctuation analysis can be obtained from his- 
tograms on x of varying bin sizes. The variance difference 
between object and reference distributions plotted as a 
function of increasing bin size exhibits a rapid initial rise 



when the bin size is less than the gaussian width, followed 
by asymptotic approach to a slower linear increase that 
continues indefinitely. That fluctuation scale dependence 
is the integral on the difference variable of the gaussian 
in the autocorrelation plus the uniform background, and 
may be processed by the methods described in this paper 
to infer the average properties of the parent distribution, 
including the gaussians. 

E. Paper outline 

We first derive the exact linear relationship between 
the bin-size dependence of fluctuations measured in a 
single-particle space and the integral of an autocorrela- 
tion projection of the two-particle difference histogram. 
We then introduce an inversion method to solve the in- 
tegral equation so that given the fluctuation information 
the autocorrelation may be inferred. Finally, we demon- 
strate physical interpretation of autocorrelations and the 
precision of the inversion method using Monte Carlo sim- 
ulations relevant to nuclear collisions. 



IV. SUMMARY OF THE METHOD 

We introduce the basic concepts of the paper and de- 
fine the integral equation between fluctuation scale de- 
pendence and two-particle correlations in the form of au- 
tocorrelation distributions. 

Spaces and measures The particle momentum in a 
relativistic nuclear collision is expressed as p — > (ptiVi <^)j 
where pt is the transverse momentum relative to the col- 
lision axis z, 77 = 5 ln{ ~^ } is the pseudorapidity and </> 
is the momentum azimuth angle. Given two-particle mo- 
mentum space (p\,p2), each momentum component x, 
bounded by some detector acceptance Ax, defines a two- 
particle subspace (x\,X2). To accommodate numerical 
integration we define a two-tiered binning system includ- 
ing microbins and macrobins. Subspace x is partitioned 
into microbins of fixed size e x (limited below by detector 
resolution). Macrobins, with variable bin size or scale Sx 
such that e x < Sx < Ax, are defined as combinations of 
the underlying microbins. Variation of fluctuations with 
macrobin scale is the basis of the fluctuation inversion 
analysis. 

If a measure W is distributed on x the corresponding 
histogram bin contents are denoted w (e.g., scalar sum 
of particle pt or particle multiplicity n in a bin) . We refer 
to single-particle measure density p(w; x) and histogram 
bin contents w a ~ e x p(w; x a ), or two-particle density 
p(w; xi, X2) and histogram w a Wl ~ e x p(w;x a ,Xb). An 
object histogram is formed with sibling pairs taken from 
same events. The reference histogram is formed with 
mixed pairs taken from different events. Given object 
distribution p bj and reference distribution p re f the dif- 
ference density (correlation function) Ap = p bj — p re f 
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represents the desired correlation structure, and Aw = 
w bj — w re f represents the corresponding difference his- 
togram. 

Fluctuations In a nuclear collision a total weight 
W (e.g., total transverse momentum P t or multiplicity 
N) falls within a detector acceptance (At?, A<p) divided 
into M macrobins of size (Sn,84>). Each macrobin (rep- 
resented by single scale Sx for brevity) contains in each 
event some integrated weight w(6x) as a random vari- 
able. The conventional per-bin variance is o~^(8x) = 
{w(Sx) — w(Sx)} 2 , where the overlines represent aver- 
ages over all macrobins of size Sx in all events of an 
ensemble 0. The corresponding scale-dependent per- 
particle variance is defined by a^^Sx) = a 2 v (Sx)/fi(Sx). 
Its small-scale limit o 2 ^ (n — 1) is ensemble- mean single- 
particle pt variance <r| t if w — p t , and is defined to be 
1 (Poisson limit) if w = n. The scale-dependent vari- 
ance difference defined by Aa^^Sx) = a^^(Sx) — a\ is a 
definite integral of two-particle correlations over a scale 
interval ^(j; as demonstrated below. 



Autocorrelations For distribution f(x) the auto- 
correlation is A(x&) = f(x) ■ f(x + xa), where the over- 
line denotes an average over some bounded interval (ac- 
ceptance) on x. For relativistic nuclear collisions the 
2D joint autocorrelation on (77a,<M) is a nearly lossless 
projection from the 4D (771, 772, </>i, $2) pair-momentum 
space. We distinguish between autocorrelation his- 
tograms A and autocorrelation densities p(x&)- We de- 
fine histogram difference AA = A bj — A re f and corre- 
sponding density difference Ap = p bj — p re j . 

Integral equation Variation of Aer^ on bin scales 
(Sr],S(f>) is represented by an integral equation which can 
be inverted to obtain an autocorrelation on difference 
variables (?7a,0a)- The integral equation is 

Aal l/ (me rn ne (j ,) = 2a l2 Aa w /(me v ,ne <t> ) (1) 



4 2J K 



mn\kl 



Ap(w;ke v ,le t j ) ) 



which also defines difference factor Ao~ w /. Derivation of 
that integral equation and the techniques required for its 
inversion are the main subjects of this paper. 



V. FLUCTUATION SCALE DEPENDENCE 



where s is a macrobin index, M(Ax, Sx) is the event-wise 
number of occupied macrobins of size Sx in acceptance 
Aa;, and the bar denotes an average over all events. This 
quantity is also a bin-based estimate of the correlation 
integral C2(w; Ax, Sx) [T^ . Because the total variance is 
extensive with respect to detector acceptance Ax we de- 
fine several intensive forms which result in, respectively, 
per-bin and per-particle variances and a total-variance 
density on x 



Y? w {Ax,5x)/M(Ax,8x) = a 2 w {Sx) 



(3) 



Y>i{Ax,5x)/N{Ax) = {Aw) 2 /n 



El(Ax,Sx)/Ax 



dn/dxa^^Sx) 



Those factorized forms separate acceptance Aa: and scale 
Sx dependencies. The optimum choice among them de- 
pends on the physical source(s) of nontrivial fluctuations. 

Variance differences In the small-scale limit where 
the particle occupancy of all occupied bins is 1 we have 
S^(Ax, Sx <§; Ax) — > N(Ax) cr|, the inclusive single- 
particle variance of w. For central limit conditions (in- 
dependent w samples from a fixed parent distribution) 
the total variance is scale invariant: Y, 2 a (Ax,Sx) = 
N(Ax) (7^ over any scale interval for which CLT con- 
ditions are valid. The scale invariance of total variance is 
an alternative statement of the central limit theorem |10| . 
Deviations of S^(Ax, Sx) from CLT invariance indicate 
the presence of w correlations. That connection leads to 
the definition of a fluctuation measure which is differen- 
tial across a scale interval: the total-variance difference 

AZ 2 w (S Xl ,Sx 2 ) ee Zl(6x 2 )-Zl(8 Xl ), (4) 

which is zero if CLT conditions apply over the scale in- 
terval. The total variance difference for 5x 2 = Sx and 
5x\ <C Aa; (n — 1) is given by 



AS; 



,(Ax,Sx) = N(Ax) {al/{Sx) - al) 



N(Ax)Aa 2 w/ (5x), 



(5) 



which defines per-particle variance difference Aa^i. The 
variance difference across a scale interval is a measure of 
correlations in p(w; x) integrated over that scale interval; 
it is equivalent to the integral of a difference autocorre- 
lation across the same scale interval, as described below. 
From that equivalence we obtain the integral equation re- 
lating autocorrelations to fluctuation scale dependence. 



Variances We define the scale-dependent total vari- 
ance of difference histogram Aw(Sx) by 



M(Ax,6x) 

ZI(Ax,5x)ee {^s(Sx)} 2 , (2) 

s=l 



VI. AUTOCORRELATIONS 

The autocorrelation of a two-particle distribution on 
binned space (xi,X2) consists of averages over diagonal 
regions of the 2D space which then occupy single bins of 
the autocorrelation histogram on the difference variable. 
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We describe the construction of autocorrelations from a 
two-particle distribution by two direct methods (i.e., by 
pair counting, not fluctuation inversion). 

Two-particle spaces In Fig.^ two-particle density 
p(w; x%, x%) is distributed on space (x\,X2). It is conve- 
nient to define sum and difference variables xs = x\ +X2 
and Xa = Xi~X2- Distance along those variables is mea- 
sured on x as shown in Fig. Q] to avoid factors \/2. If the 
correlation structure of single-particle density p(w; x) is 
independent of mean position on x or stationary (equiv- 
alently, p(w; x\ , X2) — > p(w; xs, xa) is independent of 
sum variable xs) the corresponding autocorrelation den- 
sity on xa represents a lossless projection (all struc- 
ture retained) of the two-particle density onto a lower- 
dimensional space: p(w; x^, xa) — * p(vj;xa)- For rel- 
ativistic nuclear collisions that symmetry is nominally 
exact on cj) (assuming ideal detection) and approximate 
on 77 (depending on the validity of 'boost invariance' over 
some pseudorapidity interval). 




FIG. 1: Integration schemes for direct construction of auto- 
correlations: 1) normalized strip integrals (left panel) and 2) 
diagonal bin averages (right panel). Macrobins have variable 
size Sx, microbins have fixed size e x . Bins on difference vari- 
able xa are indexed by integer k. 



Autocorrelation histograms Autocorrelation his- 
tograms can be constructed by three methods: 1) more 
directly (Fig. ^ left panel) by sorting particle pairs into 
ID bins on momentum difference variables xa as nor- 
malized strip integrals Q , 2) less directly (Fig. ^ right 
panel) as weighted averages of diagonal combinations of 
microbins from a 2D histogram on space (x\,X2), 3) in- 
directly by inverting the scale dependence of a moment, 
variance or variance difference. 

We introduce for the numerical integration microbins 
of fixed size e x such that there are m$ microbins in a mac- 
robin of variable size (e x < Sx < Ax) and toa microbins 
in a fixed acceptance of size Ax. If difference variable xa 
is partitioned into microbins indexed by k, two autocor- 
relation types, which depend on binning strategy on the 
difference variable, can be defined for an aperiodic distri- 
bution on x. Those autocorrelation types correspond to 
the two direct construction methods described above and 
illustrated in Fig. Q] (left and right panels respectively): 



1) square bin edge on zero (strip integral - particle pairs 
indexed by i,j falling within the k th xa bin), then 



1 N 

A * {W) S m*-\k\ + l/2^ {WiWi)iiek > (6) 

ij=i 

with |fc| <E [1,toa] {wi,Wj are single-particle values); and 
2) triangular bin centered on zero (averages of 2D mi- 
crobin contents indexed by a, 6), then 

^ ma- |fe| 

4 2) H = ^3|fc| E ( 7 ) 

1 1 a— 1 

with |fc| £ [0, 777a — 1] (average over all bins on diagonal 
bin array within acceptance, Fig. ^ (right panel) with 
Sx — » Ax). Overlines indicate averages over all events. 
For a distribution periodic on x the type-1 autocorrela- 
tion is given by A^\w) = l/rriAY^=i( w i w j)i3ek, with 
particles of any pair i,j falling within one period on x 
(the period then defines Ax). Binning strategies and in- 
dexing on xa for the two cases are shown in Fig. [21 For 
autocorrelations constructed from a particle list (most 
efficient) the square bin index applies. For prebinncd 
2D space (xi,x%) the triangular bin index applies. For 
fluctuation inversion leading to autocorrelations either 
method may be used, but the rectangular bin index fa- 
cilitates direct comparison between fluctuation inversion 
and the more efficient direct method 1) of autocorrelation 
construction. 

k = 1 2 rectangular bin index m 

k = 1 2 triangular bin index m-1 8x X A 

FIG. 2: Binning schemes on difference variable xa- 1) square 
bins on xa (Fig. left panel) with |fc| £ [l,m], and 2) 
triangular projections of 2D bins (Fig. Q right panel) with 
\k\ £ [0,m - lj. 



Difference autocorrelations AAk — Ak.obj — 
Ak,ref represents a subset of the correlation content of 
an object distribution. The type-1 number 77 reference 
autocorrelation obtained with mixed pairs (particles se- 
lected from different but similar events) is A^\ e An) = 

1/(toa - \k\ + l/2)X)ijeft = n l (hJ are particle in- 
dices). The type-2 equivalent based on bin averages is 

4? e /W = 1/Va- |fc|)E:=r' fc| ^"a +fc = n\. The 
type-1 difference autocorrelation for measure w is de- 
fined in terms of sibling (same-event) and mixed pairs 
as AAj^H = l/(m A - \k\ + l/2)E„ efc A^A^-. That 
form is illustrated in Fig. [T] (left panel). We can utilize bin 
averages w a ,ref = w a to define type-2 difference autocor- 

relation AA%\w) = l/(m A - I^DE^i"'* 1 Aw * Aw a+k . 
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(2) 

As defined, AA k (w) is a projection by averaging onto 
its difference variable of binned two-particle covariance 
distribution Aw a Aw b on (a; 1,2:2) i a i b are microbin in- 
dices). That form is illustrated in Fig. ^ (right panel, 
with Sx — ► Ax). The difference autocorrelation, being 
a projection of A.p(x\,X2) onto its difference variable, is 
also a covariance distribution (most obvious for the type- 
2 formulation above) , except for the k = bin which is a 
variance at the microbin scale. Figs. (right panels) and 
[5] show examples of difference autocorrelations obtained 
from fluctuation inversion. 



VII. THE INTEGRAL EQUATION 

We next obtain an exact correspondence between the 
per-particle variance and a type-2 difference autocorre- 
lation. We use the relation between binning schemes in 
Fig-Glto express per-particle variance as the scale integral 
of a type-1 difference autocorrelation, most efficiently ob- 
tained from a direct autocorrelation analysis. We then 
generalize to 2D momentum subspace (77, <ft) to obtain 
the required integral equation relating fluctuation scale 
dependence on (Srj, S(f>) to joint autocorrelations on axial 
difference variables (j]a, <M)- 

Variances and type-2 autocorrelations Starting 
with the definition of total variance we obtain the exact 
relation of the per-particle variance to a weighted sum 
of type-2 autocorrelation elements. The total variance 
as a sum over occupied macrobins is given in Eq. 
Since Aw s (Sx) — J2 a£s A.w a (e x ), the corresponding per- 
particle variance can be written as a sum over 2D mi- 
crobins 



m$—l b—a+k 



a -/ {5x) = Wx) ^ ^ ^ a Aw b , (8) 

fe=l— mg a,b—l---rns 

where the sums over indices a, b include only microbins 
within macrobins at scale Sx, as in Fig. ^ (right panel), 
and the upper overline indicates an average over all oc- 
cupied macrobins in the acceptance. Rearranging and 
noting that n(8x) = m$n(e x ) we obtain 



°w/(& x = m * £ z) 



m s -l 1 , 1 

ms - \k\ 



k—l — ms 



Ills 



(9) 



n k (e x ) 
n(e x ) 



ms 



\k\ 



b—a-\-k 

E 



\Jn a n h ( Aw a Aw b 



where n k (e x ) = 1/ (ms - |fc|) 2 0) f, a+ V n a n b is a geomet- 
ric mean over microbins within macrobins at scale Sx, 
and the upper overline denotes an average over all mac- 
robins in the acceptance. 

The quantities in parentheses are two-particle 
difference-histogram elements equivalent to density ra- 



tio e x {Ap(w) / p re f(n)}ab- That 2D histogram rep- 
resent per particle correlations in measure w on space 
(xi,X2). The quantity in square brackets is a weighted 
average of those histogram elements over the k th diag- 
onal which estimates type-2 difference autocorrelation 



AA^\w)/ \J A y £ re f(ri), the per-particle autocorrelation 
distribution on xa- The quantity in curly brackets es- 



i(2) 



timates autocorrelation elements AA^ (w) / \J A y k ' re f (n) 
which would be inferred from inversion of fluctuation 
scale dependence (cf Eq. JTTJ). ^ quantity n k varies sig- 
nificantly (possible on rj A , depending on collision system, 
detector acceptance and required precision) the 0(1) 
multiplicity-ratio factor n k /n may be removed as a sys- 
tematic correction to obtain the hatted autocorrelation 
in square brackets. This per-particle relation offers an 
optimized match between two-particle correlations and 
fluctuation measures. 

Variance differences and type-1 autocorrela- 
tions We now transition from type-2 to type-1 autocor- 
relations in order to relate fluctuation scale dependence 
to the more efficient direct method of obtaining autocor- 
relations. The quantity in curly brackets in Eq. esti- 
mates the typc-2 difference autocorrelation defined pre- 
viously as an average over all microbins on the k th diago- 
nal within the acceptance (Fig.^ right panel, Sx — > Ax). 
We rewrite Eq. JHJ), replacing the quantity in curly brack- 
ets by the difference autocorrelation which it estimates. 
We then re-express the variance difference in terms of an 
equivalent expression for the type-1 difference autocorre- 
lation 

2 / \ "sr m s -\k\ { AA k {w) ),.,„,. 
a w /(m s e x ) = ^ — — — { /a } (10) 



,(2) 



k—l — mx 



m s I yJA k , ref (n) 



2 gm 5 -fc + l/2 J AA k (w) 

k=l 



(1) 



^/A k:ref (n) 



m s - fc + 1/2 x 

k=l 



m s 



n k (e x ) 
n(e x ) 



1 



N 



1/2 E 



Aw, Ai 



m A - k + 1/2 n k (e x ) 



where Awi = wt — w for particle i, w is the inclusive mean 

of w, nl(e x ) = Y.i!jek/( m ^ ~ k + 1/2) for mixed pairs, 
and the quantity in square bracket is the autocorrelation 
obtained by direct method (1). Either expression (1) or 
(2) above can be used to define the integral equation re- 
lating fluctuation scale dependence to difference autocor- 
relations. The two forms correspond to the two binning 
schemes on the difference variable shown in Fig. |2j 

Inclusive variance <7^ in variance difference Aa^^Sx) 
corresponds to self pairs in the pair sums of Eq. l|lUfl • If 
those pairs are excluded from two-particle pair distribu- 
tions and the autocorrelation definitions then the rela- 
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tions above apply to the variance difference, rather than 
the scale-dependent variance. For the type-1 autocorre- 
lation we then obtain the integral equation 



Ao i/(w£x) 



2^ K m -i 

k=l 



AA k [w] e x ) 
' \/A k ^ e f(n; e x ) ' 



(11) 



where K m: k = (to — k + l/2)/m is an integration kernel 
representing the macrobin scheme, and we have dropped 
the subscript S from index to. Autocorrelation difference 
density Ap is related to the autocorrelation difference 
histogram by AAk(w;e x ) = e 2 Ap(w; k e x ). Likewise, 
A-k,ref(n]e x ) = e 2 p re f(n; k e x ). The integral equation 
is then re-expressed in terms of a density ratio as 



fc=i 



rn ; k 



Ap(w; k e x ) 
y/pref(n;ke x ) 



(12) 



2D integral equation For a 2D scaling analysis on 
(rj,(f) we generalize 5x — > (5r],S4>) to obtain the per- 
particle w variance difference (also defining difference fac- 
tor Aa w / ,11]) as a 2D discrete integral equation 

AtT ™/( TOe »?' ne$) = 2(jy,Aa w /{me n ,ne (t> ) (13) 



4 ^ e^£0 K mn -ki 

k,l=l 



Ap(w;ke n ,le <t> ) 
y/pref(n;ke v ,le l j > y 



with kernel K mn -ki = (m — k + 1/2)/to ■ (n — I + l/2)/n 
representing the 2D macrobin system. This is a Fredholm 
integral equation which can be inverted to obtain auto- 
correlation density ratio Apu / ^Jp re f,ki as a per-particle 
correlation measure on difference variables (to\, 4>a) from 
fluctuation scale dependence of Aa^^Sr/, 5<j)). 



VIII. THE INVERSE PROBLEM 

We have established the relationship between a vari- 
ance difference and the integral of an autocorrelation. 
Given the former as experimental data we wish to in- 
fer the latter as conveying the underlying physics in a 
more interpretable form. Joint (2D) autocorrelation dis- 
tributions on difference variables (^Aj^a) are inferred 
from fluctuation scale dependence on (chy, Sip) by invert- 
ing those distributions according to Eq. (|13fl (those au- 
tocorrelation distributions are by construction symmetric 
about the origin on (t^a, <M))- I n standard nomenclature 
the variance difference on the left in Eq. (|13fl is 'data' D 
and the autocorrelation density ratio on the right is 'im- 
age' I. Schematically, D = TI <-> Acr^ = KAp/^/p. 

Extracting the image Inverting an integral equa- 
tion from data to image is usually nontrivial and has en- 
gendered a large literature |12( . Integration and inversion 
are mappings (matrix transformations T, T _1 ) from one 
vector space to another. Forward integration from image 



to data can be represented by D = T I + N, (N repre- 
sents noise added to the data as part of measurement) 
with formal inverse I' = T D = I + T _1 N, where I is 
the true image. The image produced by naive mapping 
I' = T 1 D is typically dominated by small- wavelength 
components of image noise T _1 N, since T _1 is a form of 
differentiation (acting as a 'high-pass' filter). A smooth- 
ing procedure called regularization is used to reduce im- 
age noise while minimizing distortion of the true image. 



Regularization By smoothing the inverted image 
I' a = T" 1 D with a regularization procedure (smoothing 
parameter a is defined below) we control the impact on 
the image of the statistical noise on the data. We treat 
each bin of the image as a free fitting parameter, but add 
a smoothing term with Lagrange multiplier a to the fit 
X 2 to obtain xl = ||D - TI^|| 2 + a||LI^|| 2 . The first 
term measures the deviation of the integrated image from 
the data. The second term, with linear operator L, mea- 
sures the small- wavelength A structure of the image. The 
optimum choice of a maximally reduces small- A noise on 
l' a without significantly distorting the signal - the mean- 
ingful structure on true image distribution I. Operator L 
can take various forms The underlying assumption 
of a regularization procedure is that the statistical 'noise' 
is sufficiently distinct from the true correlation 'signal' 
on wavelength A that operator L, a measure of small-A 
structure on 1^, can distinguish the two. We eliminate 
most of the small-A noise while retaining as much large-A 
signal as possible. As a increases, more small-A noise is 
eliminated. At some point increasing a reduces the sig- 
nal. By comparing the two elements of xL as m Fi§- El 
(left panel), the optimum value of a can be determined. 

Statistical and systematic errors The autocor- 
relation (image) error distribution has two components: 
statistical error which survives smoothing and system- 
atic error due to image distortion by smoothing. The 
optimum value of a should adequately reduce statisti- 
cal noise relative to signal while minimally distorting the 
signal. The residual statistical noise and the signal dis- 
tortion then depend on the optimum value of a. The 
linear relation between statistical errors on data and im- 
age is similar to that for data and image themselves. It 
is straightforward to invert that relation (using the op- 
timized a value) to obtain mean-square errors S 2 I a for 
the image histogram. Because of the structure of differ- 
ence autocorrelation density ratio Ap/ ^Jp re f those errors 
tend to be uniform on the difference- variable space, mak- 
ing their characterization simple. 

To estimate smoothing distortions, inferred image l' a 
can be integrated forward to become a new data distribu- 
tion, and then inverted again to obtain a twice-smoothed 
image. That procedure is represented schematically by 

1) l' a = T^D (smoothed image from data inversion), 

2) D„ = TI' Q (forward-integrated smoothed image), 3) 
l' aa = T t ^ 1 D Q (twice-smoothed image). By compar- 
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ing image distributions 1^ and V aa , assuming that the 
smoothing process is linear, one can estimate the smooth- 
ing distortions, usually largest at points of maximum gra- 
dient on the image. The smoothing error can then be 
displayed as a histogram or summarized as a fractional 
error of peaked structure on the image distribution. 



IX. MONTE CARLO EXAMPLES 

We apply the inversion method to Monte Carlo data for 
three reasons: 1) these Monte Carlos produce event en- 
sembles with fairly simple correlation content which are 
amenable to the type of fluctuation/correlation analysis 
we wish to apply to nuclear collision data, 2) inversion 
of Hijing heavy ion fluctuations illustrates how the in- 
version process leads to simple and direct interpretation 
of fluctuations in terms of physical mechanisms and 3) 
comparison of the Pythia fluctuation inversion with a di- 
rect computation of Pythia correlations (easily done for 
p-p collisions) confirms that the inversion process pre- 
cisely reveals the underlying autocorrelation with negli- 
gible distortion. Those basic demonstration goals could 
be accomplished as well with a simpler physics Monte 
Carlo, or with a Monte Carlo having no physical cor- 
respondence to data {e.g., the gaussian on a flat back- 
ground). The point of this exercise is to make controlled 
comparisons between inversion and direct pair-counting 
analysis to establish the precision of the fluctuation in- 
version, and to perform a controlled inversion study of a 
Monte Carlo with known physics content for calibration, 
before applying the method to real data. 

Hijing Au-Au 200 GeV p t correlations Event 
generator Hijing-1.37 [lg was used to produce a total 
of 1M minimum-bias Monte Carlo collision events of two 
types: 1) quench-off Hijing - jet production enabled but 
no jet quenching and 2) default or quench-on Hijing - jet 
production and jet quenching both enabled. Each event 
type was then split into six centrality classes consisting 
of roughly equal fractions of total cross section. A more 
complete analysis of Hijing fluctuation data and inversion 
to pt autocorrelations is described in [!j • For the fluctua- 
tion analysis charged particles with pseudorapidity \r/\ < 
1, transverse momentum p t <E [0.15, 2] GeV/c and full az- 
imuth were accepted. Fluctuation scale dependence was 
then determined using the variance difference defined in 
Eq. JSJ, with w — > (p t — npt) and w / — > pt : n. 

Fig. El shows fluctuation scale dependence and cor- 
responding autocorrelations for 0-15% central collisions 
and two Hijing configurations: quench-off (top pan- 
els) and quench-on (bottom panels). In the left panels 
are variance differences Aer 2 t . n (J?7, <5</>) [units (GeV/c) 2 ], 
which typically increase monotonically with pseudorapid- 
ity scale (5r)) but are more complex on azimuth scale 
(S(f>). Autocorrelation distributions Ap/ y/p re f [also units 
(GeV/c) 2 ] on (t7a,</>a) shown in Fig. (right panels) 
contain two features: a same-side component (|</>a| < 




FIG. 3: Left panels: Act 



2 

pt-n 



(GeV/c) 2 distributions on bin- 
ning scales (Sn, 8(f>) for central 0-15% collisions and two Hijing 
states: quench off (top panel) and quench on (bottom panel). 
Right panels: Corresponding autocorrelations on difference 
variables (»?a,0a). 



7r/2) and an away-side component (7r/2 < \4>/\\ < tt)- 
The dominant same-side peak can be identified with 
minimum-bias parton fragments (minijets): no high-p f 
trigger condition is imposed in the analysis |9(. 
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4: Left panel: Distribution of ||D — TI a || 2 (dots) and 
| 2 (triangles) vs a for the 0-15% central quench-off data 
in Fig. |21 (upper-right panel). Deviation of ||D — TI a || 2 from 
the extrapolation (horizontal solid line) indicates systematic 
error due to image smoothing. The vertical dashed lines 
bracket the optimal choice of a. Right panel: Distribution 
of estimated smoothing error Iq, — I aa on (j7a>0a), also for 
the 0-15% central quench-off autocorrelation in Fig. |21 (upper- 
right panel). Note the 10 x scale reduction compared to Fig. [31 



An illustration of a optimization is given in Fig. 01 
Distributions ||D-TI Q || 2 and ||LI Q || 2 on a shown in 
Fig. 0] (left panel) can be used to determine the optimum 
a value. The ||D — TI Q || 2 distribution is nearly con- 
stant at smaller a, indicating no loss of valid structure 
by smoothing. The ||LI Q || 2 distribution falls rapidly, in- 
dicating substantial noise reduction. At some value of a 
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||D — TI Q || 2 begins to rise from the constant (extrapo- 
lated by the solid line), the deviation representing infor- 
mation loss from the image. Beyond that point there is 
a tradeoff between smoothing distortion and noise reduc- 
tion. In the example of Fig. 0] the noise exhibits a sta- 
tionary point at the same position (upper dotted line). 
The optimum choice of a is in the region bracketed by 
the vertical dashed lines. 

Autocorrelation errors have two components: statisti- 
cal errors derived from corresponding errors in fluctua- 
tion data, and distortion errors due to smoothing. Statis- 
tical errors on the image are found by inverting statistical 
errors on fluctuation data. The structure of this error dis- 
tribution is usually simple, and can be characterized by 
a single number for each autocorrelation. Smoothing er- 
rors are estimated by passing inferred image I a a second 
time through the integration/inversion process to obtain 
l aa - Difference l a — I aa provides an upper limit on the 
smoothing distortion (some of the change includes further 
reduction of statistical error) . A typical image-smoothing 
error distribution is shown in Fig. 0] (right panel). The 
per-bin systematic error on autocorrelation distributions 
is typically a few percent of maximum autocorrelation 
values, somewhat larger peak errors lying within small 
regions corresponding to large image gradients. 

These Hijing results illustrate the process of reduc- 
ing event-wise fluctuations to autocorrelations and their 
physical interpretation. The autocorrelations revealed in 
this study and in [9j have features consistent with the 
physical model of jet production and quenching in that 
Monte Carlo. They represents the first determination of 
transverse momentum correlations from (p t ) fluctuations 
for any nuclear collision data or model. Transverse mo- 
mentum autocorrelations provide a much more critical 
and differential test of the Hijing model than was previ- 
ously conducted (the model fails those comparisons with 
data). The inversion process clearly reveals what mech- 
anism in Hijing generates (pt) fluctuations, which was 
previously a matter of speculation. These jet correla- 
tions were not obtained with a jet trigger per se (leading- 
particle or jet-cone analysis), but instead result from an 
unbiased, symmetric two-particle correlation analysis. 

Pythia p-p 200 GeV n correlations The Pythia- 
6.131 Monte Carlo 0] was used to provide a high- 
statistics sample of 200 GeV p-p collision events. Anal- 
ysis of those data provides a demonstration of the pre- 
cision possible with fluctuation inversion and a check on 
the algebraic system of this paper. We make a compari- 
son between autocorrelations determined directly accord- 
ing to Eq. JBJ and indirectly by fluctuation inversion ac- 
cording to Eq. Ql.'jfl . In Fig. are charge-independent 
(left panel, like-sign + unlike-sign pair combinations) 
and charge-dependent (right panel, like-sign — unlikc- 
sign pair combinations) pair-number autocorrelations on 
difference variables (?/a,<M) 0- The autocorrelations 
were obtained by fluctuation inversion (upper panels) 
and directly from binned pair numbers (lower panels). 



Detailed comparison of the shapes of each pair of dis- 
tributions indicates that fluctuation inversion is a very 
precise method of determining autocorrelations. 




FIG. 5: Pair- number autocorrelations on difference variables 
(f?A,<M) f° r P"P collisions at 200 GeV from Pythia-6.131 ob- 
tained from fluctuation inversion (upper panels) and from di- 
rect pair counting (lower panels) for isoscalar (left panels) and 
isovector (right panels) autocorrelations. 



The shift of vertical scale for upper and lower panels in 
Fig.|S]is due to the choice of sibling vs mixed pair normal- 
ization used in the direct determination of pair-number 
autocorrelations. Pair-number autocorrelations are con- 
structed from pair sums hij = n^j JV ■ m$ normalized to 
unit integral for sibling and mixed pairs [3j . Normalized 
pair-ratio histograms fy = n.ij,abj/nij,ref ar e first formed 
on microbinned spaces (771,772) and (4>\,4>2)- Those dis- 
tributions are then projected to form type-2 autocorrela- 

tion ratios AA y k /Ak : ref = ^fe — 1 on (rjA, 4>a)- The pair- 

(2) 

number difference autocorrelation AA k so constructed 
has zero integral, and therefore always varies about zero 
as in the lower panels of Fig. [5J Autocorrelations from 
fluctuation inversion retain their absolute vertical posi- 
tion: the offset missing from the zero-integral type-2 au- 
tocorrelation is equivalent to a fluctuation measurement 
at the acceptance scale. The factor N in Fig. [5] defines a 
transitional form: NAA k ' /Ak t ref — AnA<f) ■ Ap/ ^/p r ef 

The fluctuation inversion procedure is a general anal- 
ysis technique. It is not unique to these Monte Carlos or 
to nuclear physics. It can be applied to any problem in 
which correlated distributions are produced as an event 
sequence and there is persistent event-wise structure over 
an event ensemble. 
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X. SUMMARY 

Event-wise fluctuations in the hadronic final state of 
in nuclear collisions, a manifestation of multiparticle cor- 
relations produced by collision dynamics, are relatively 
easy to calculate but often difficult to interpret. Auto- 
correlations provide a powerful method for direct study 
of multiparticle correlations but are computationally ex- 
pensive when constructed directly from pair ratios. For 
relativistic collisions, autocorrelation projections to sub- 
spaces of reduced dimensionality retain most of the cor- 
relation information of the full two-particle momentum 
space in a more compact form. 

In this paper we have derived the exact relation be- 
tween autocorrelations and fluctuation scale dependence 
in the form of an integral equation. We have introduced 



the method required to invert that equation, and the 
procedure used to solve the technical problem of regular- 
izing statistical noise. Using Monte Carlo simulations we 
have made detailed comparisons of autocorrelations de- 
termined directly from pair ratios and from fluctuation 
inversion with regularization from the same particle data 
to illustrate the precision of the method. 

The analysis of the Hijing Monte Carlo illustrates how 
inversion of fluctuations to autocorrelations permits di- 
rect interpretation in terms of physical mechanisms, in 
this case angular correlations of transverse momentum 
from jet fragmentation. The inversion method applied 
to transverse momentum fluctuations provides the first 
access to velocity correlations in nuclear collisions. This 
procedure is not unique to nuclear physics data, but can 
be applied to any point sets arranged in event ensembles. 
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